Four-photon Entanglement as Stochastic-signal Correlation 
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A fully classical model of a recent experiment exhibiting what is interpreted as teleportation and 
four-photon entanglement is described. It is shown that the reason that a classical model is possible, 
contrary to the current belief, results ultimately from a misguided modification by Bohm of the EPR 
Gedanken experiment. Finally, teleportation is reinterpreted as the passive filtration of correlated 
but stochastic events in stead of the active transfer of either material or information. 

PACS numbers: 03.65.Bz, 03.67.-a, 41.10.Hv, 42.50.Ar 



O 
O 

> 
O 

(N 
(N 

> 

m 
O 
in 
o 



^— > ■ 



In a recent letter J.-W. Pan et al. described a 
demonstration of the experimental observation "of pure 
four-photon GHZ entanglement produced by paramet- 
ric down-conversion and a projective measurement." jlj 
They add that this experiment "demonstrates teleporta- 
tion with very high purity" and that "the high visibility 
not only enables various novel tests of quantum nonlocal- 
ity, [but] it also opens the possibility to experimentally 
investigate quantum computation and communications 
schemes with linear optics." This demonstration was 
achieved using a novel and ingeniously simple (in con- 
cept, not necessarily in concrete realization) setup em- 
ploying two sources of entangled photons which feed the 
two faces of a polarizing beam splitter (PBS). The au- 
thors further claim that this setup gives results that fi- 
nally preclude all doubt that nonlocal effects ensue from 
quantum entanglement. 

It is the purpose of this letter to contest this last claim. 
This will be achieved by proposing a classical model of 
their experiment, without 'nonlocality,' which fully du- 
plicates its results. 

First, we briefly review the experimental setup. (See 
Fig. 0) Two independent entangled photon pairs are cre- 
ated by down-conversion in a crystal pumped by a pulsed 
laser. The laser pulse passes through the crystal creating 
one pair (A), then is reflected off a movable mirror to 
repass through the crystal in the opposite direction cre- 
ating a second pair (B). One photon from each pair is fed 
directly through polarizers to photodetectors (photons 1 
and 4). The other photons (2 and 3) are directed to op- 
posite faces of a PBS, (i.e., a beam splitter which reflects 
vertically and transmits horizontally polarized photons) 
after which the exiting photons are sent through vari- 
able polarizers into photodetectors. The path lengths for 
photons 2 and 3 are adjusted so as to compensate for the 
time delay in the creation of the pairs. By moving the 
mirror, the compensation can be negated to permit ob- 
serving the disappearance of interference caused by lack 
of simultaneous "cross-talk" between channels 2 and 3. 

The reported observations are the following: Of all the 
16 possible polarizer settings regimes for which 9 n = 
or ir/2 w.r.t. the horizontal axis of the PBS, only 
{0, tt/2, tt/2, 0} and {tt/2, 0, 0, tt/2} yield a (sub- 
stantial) four-fold coincidence count: C; the regime 
{7r/4, 7t/4, 7t/4, 7r/4} occurs with an intensity C/4.2 and 




FIG. 1: Schematic of the experimental setup for the mea- 
surement of four-photon GHZ correlations. A pulse of laser 
light passes a nonlinear crystal twice to produce two entangled 
photon pairs via parametric down conversion. Coincidences 
between all four detectors are used to study the nature of 
entanglement. 



the regime {7r/4, 7t/4, 7r/4, — 7r/4} with zero intensity. 
Further, both of the latter regimes yield an intensity of 
C/8 when photons 2 and 3 do not overlap. 

The model proposed herein, like all others, is based on 
certain assumptions, which, being quite different from 
those couched in the notation and vocabulary of QM, 
must be delineated explicitly. They consist of the follow- 
ing: 

1. Electromagnetic radiation is comprised of continu- 
ous waves as described by Maxwell's equations. 

2. All detectors in the optical region of the electro- 
magnetic spectrum exploit the photoelectric effect. 
They convert continuous radiation to an electron 
current. Electrons are discrete objects; their gen- 
eration in a photodetector effectively digitizes the 
signal associated with incoming radiation, thereby 
evoking the impression that the radiation was itself 
somehow digitized into units (photons). The last 
inference is unwarranted; we can not know what 
form the incoming energy actually had, we are re- 
stricted to inferring its nature from the electron 
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current. It is known empirically that photoelec- 
trons are ejected from a photodetector randomly 
but in proportion to the energy density, E 2 , of the 
incoming signal. For the model, we take the con- 
version efficiency to be 100%. 

The nonlinear crystal generating signals by para- 
metric down-conversion is taken to be a "black box" 
which, by virtue of its structure emits randomly 
pairs of signal pulses confined to the pure verti- 
cal and horizontal polarization modes. The pairs 
are anticorrelated with respect to polarization and 
each member of a pair is directed (in the well known 
manner) into two intersecting cones; all horizontal 
emissions are into one cone, all vertical into the 
other. The two variants of pairs are taken to be 
equally likely. "Entangled" radiation samples ex- 
tracted from the two points where the cones in- 
tersect, therefore, constitute a mirror image pair 
of random sequences of individual electromagnetic 
pulses of both pure modes. 

The intensity of four-fold coincidence detections 
among four photodetectors is calculated using non 
quantum coherence theory. Coincident count 
probabilities, for a system with N, (herein N = 4 ), 
monitored exit ports are proportional to the single 
time, multiple location second order (in intensity) 
cross correlation, i.e.: 
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< n»=i E*(r n ,9 n ) ULn E(r n ,6 n ) > 
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(1) 



For this model the denominator consists of constants of 
factors of the form (cos 2 (a) + sin 2 (a)) so that it equals 1. 

These assumptions, all fully compatible with classi- 
cal physics, account for all of the reported observations. 
Eq. ([j]) was implemented numerically as follows: The 
centers are assumed to emit double pulses in opposed 
directions which are anticorrelated and confined to the 
vertical and horizontal polarization modes; i.e. their po- 
larization vectors are: 



A! 


= ( c os(n|), 


sm(n|)), 
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A 2 


= (sin(n-), 


-cos(n~)), 
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Bi 


= (cos(m^), 


sm(m-)), 


(4) 


B-2 


= (sm(m^), 


-cos(m~)), 


(5) 



where n and m take the values and 1 with a flat ran- 
dom distribution. The polarizing beam splitter (PBS) is 
modeled using the transition matrix for a polarizer, 



P(0) 



cos 2 (0) cos(0)sin(0) 
sin(0)cos(0) sin 2 (0) 



(6) 



where = n/2 accounts for a reflection, and = a 
transmission. Thus, the field impinging on each of the 
four detectors is: 

E x = P{6 1 )A 1 , (7) 

E 2 = P(0 2 ){P(0)B 2 -P(n/2)A 2 ), (8) 

E 3 = P(9 3 )(P(0)A 2 -P(n/2)B 2 ), (9) 

E 4 = P(0 4 )B 1 . (10) 

Putting Eqs. (7-10) into Eq. (1) and carrying out the 
indicated average by summing on n and m and divid- 
ing by 4 as there are four combinations, one obtains the 
expression for the fourth order (in fields) coincidence cor- 
relation as a function of the polarizer settings. The result 
is cumbersome but easily evaluated numerically for var- 
ious regimes, i.e., various selections of 0„. Ideal results 
compatible to those observed are obtained. Slight devia- 
tions in the observed data appear to be due to noise. 
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FIG. 2: The upper curve shows the effect on the intensity 
of four-fold coincidences of skewing (rotating) all polarizers 
through a given angle starting from the state {vr/2, 0, 0, tt/2}. 
The lower curve shows the same effect when one of the polar- 
izers is rotated in the opposite direction. The middle curve 
shows the effect of either of these skewing schemes when the 
timing is such that the crossover signals do not arrive simul- 
taneously with the reflected signals. Note that the separation 
of the curves and values at n/4 coincide with the observations 
reported in Ref. This diagram differs from Fig. 4 in Ref. 
|Q] in that it shows the split of these regimes as a function 
of polarizer skew for fixed delay rather than as a function of 
delay for fixed skew. 

To model regimes for which the pulse pairs were gen- 
erated with a time difference such that the cross-over sig- 
nals in beams 2 and 3 could not interfere, the sum on the 
radiation contributions from individual centers is squared 
before averaging. This procedure recognizes the fact that 
the electric fields from distinct sources (A and B, in this 
case) must be added after squaring because they do not 
interfere but still deposit energy into a photodetector. 
Changing the order of 'squaring' and 'summing' for this 
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setup affects only the crossover signals as all others are in 
orthogonal polarization modes in any case. All this leads, 
again, to results exactly mimicking the effects reported. 
Additionally, the relative count intensity in other regimes 
can be computed easily; for an example, see Fig. H 
Note that on Fig. 2 the curve for the 'state' 

< 7r/4, 7r/4, 7r/4, 7r/4 > splits or separates from the 'state' 

< 7r/4, 7r/4, 7r/4, -7r/4 > at 6* = 7r/4. Pan et al. interpret 
this splitting to indicate: "that the state of, say, photon 
2 was teleported to photon 4..." and ". . .fully demon- 
strates the nonlocal feature of quantum teleportation." 
The mere fact that this splitting is faithfully modeled 
classically, using, e.g., coherent vice Fock states, casts 
strong doubt on such inferences. 

There is nothing essentially quantum mechanical in 
this model. It does not make use of any property pe- 
culiar to QM. In particular it does not in any way as- 
sume properties of "photons" and its results do not de- 
pend, therefore, on assumptions regarding detector ef- 
ficiency (effectively taken to be 100%). It constitutes 
a fully faithful classical technique for modeling experi- 
ments carried out in a portion of the spectrum admit- 
ting macroscopic devices and detailed time tracking of 
electromagnetic fields, thereby evading the peculiarities 
of photodetectors. (Note that EPR-B correlations have 
been so observed. ||) The same type of model for EPR-B 
experiments yields exactly the same expressions as those 
given by QM which violate the relevant Bell Inequality. 

In view of the conventional wisdom that it is impossible 
to comprehend phenomena involving quantum entangle- 
ment using non quantum physics it is natural to ask 
how it is possible then that a classical model can be suit- 
able. We believe that the genesis of this misunderstand- 
ing is to be found with D. Bohm.[^| The over-arcing issue 
ultimately motivating the experiment described in Ref. 
0, and further analyzed herein is the Einstein, Podol- 
sky and Rosen [EPR] argument of 1935 to the effect that 
QM is incomplete. J(| In that paper EPR argued that the 
position and momentum of two entangled particles, cre- 
ated auspiciously, could be specified exactly in spite of 
Heisenberg Uncertainty, by exploiting symmetry. With 
this argument, they hoped to show that Heisenberg Un- 
certainty was not something fundamentally new, but just 
ignorance. This being the case, they argued, there should 
exist a deeper theory, involving heretofore "hidden vari- 
ables" that would more precisely describe nature. The 
point of their whole consideration involved the ultimate 
nature of Heisenberg Uncertainty. 

Much later Bohm modified EPR's original Gedanken 
Experiment. He, for reasons of simplicity, transfered the 
EPR argument from phase space, (x, p) to another arena, 
one involving spin.[|| This too turned out to be exper- 
imentally unrealizable, but the algebraically isomorphic 
arena involving polarized light is entirely practical. Thus, 
nowadays, the EPR issue is discussed, analyzed and ex- 
plored in terms of polarization phenomena. However, 
note that because of Heisenberg Uncertainty, (x, p) do 
not commute, i.e., [x, p] = ifr, whereas the basis oper- 



ators of polarization space Eh, E p , where E x represents 
an electric field in the x direction, do commute. That is, 
there is no Heisenberg Uncertainty among polarization 
modes. This is a fact substantiated by QM itself; cre- 
ation and annihilation operators for different modes of 
polarization, commute. Bohm's transfer of venue moved 
the issue from one in which there is Heisenberg Uncer- 
tainty into one in which there is not! Experiments in the 
polarization arena can not, therefore, address the issue 
introduced by EPR; they, by logical necessity, leave it un- 
examined. It is simply impossible to investigate Heisen- 
berg Uncertainty where there is none. 

Because the basis operators of polarization space 
(a.k.a. "qubit space") commute, all polarization phenom- 
ena ultimately must be describable with non quantum 
principles. It is for this reason, that a classical model 
can explain this experiment. Alternately, this conclusion 
follows forthwith, albeit with mostly only formalistic au- 
thority, from the Optical Equivalence Theorem. [Q 

Bohm did not justify carefully his change of venue; he 
simply declared spin operators to be equivalent and sal- 
lied forth. This has been accepted, apparently, on the 
grounds that, like the basis operators of phase space, 
[x, p) , the spin operators (a x , a y ) also do not commute. 
While this is indeed true, the reason is not the same. 
Spin operators with discrete eigenvalues pertain to the 
direction in space defined by a magnetic field. In direc- 
tions transverse to the magnetic field, the expectation 
values are not discrete but oscillate out of phase. In the 
end, the reason spin operators do not commute would 
be that it is impossible to have more than one direction 
for a magnetic field or precession axis at a time. Non 
commutation, then, is not a manifestation of Heisenberg 
Uncertainty in this case, but of geometry, i.e., of the na- 
ture of an axis of rotation. Indeed, angular momentum 
operators do not commute in classical mechanics. Similar 
remarks pertain to the privileged direction defined by the 
"k" vector of an EM wave with respect to the transverse, 
i.e., polarization directions. 

Entanglement is often cited as the core of QM, an 
idea for which perhaps Schrodinger was the originator. || 
However, if entanglement is defined in terms of the 
non factorability of a wave function, as it most of- 
ten is, then it must be attributable to a correlation 
between subsystems. Obviously, if a wave function fac- 
tors into terms each pertaining to a separate subsys- 
tem, ipi{r 1)1^2^2), then the probability density com- 
puted from this product, iplip^ri, r2)ip2'4 ! i(j'i, T2) = 
('0i(f'i)V'i( r 'i))(V'2( r 2)'02(? , 2)), also factors, in which case 
it pertains to statistically independent, i.e., uncorrelated 
subsystems. If the wave function does not factor, then 
the probability density also will not factor and, clearly, 
the subsystems are simply not uncorrelated. Correlation 
need not imply nonlocality; hereditary correlation satis- 
fies all requirements posed by the physical situation. In 
experiments involving EPR style setups, including that 
described herein, correlation can be vested in the 'pho- 
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> 


restart: a:=l: b:=l: c:=0: d:=0: e:=l: 


> 


with ( linalg) : 


> 


Pro :=z->matrix ( [ [cos (z) "2, cos (z ) *sin (z) ] , [sin ( z) *cos(z) , sin (z) "2] 


> 


7\ 1 

Al 


=vector ( [cos (n*Pi/2) , sin (n*Pi/2) ] ) : 


> 


A2 


=vector ( [sin(n*Pi/2) , -cos (n*Pi/2) J ) : 


> 


Bl 


=vector ( [cos(m*Pi/2),sin(m*Pi/2) ] ) : 


> 


B2 


=vector ( [sin(m* Pi/2), -cos(m* Pi/2) ] ) : 


> 


El 


=multiply (Pro (zl) ,A1) : 


> 


E2 


=multiply (Pro (z2) , -multiply (Pro (Pi/2) , A2 ) +multiply (Pro (0) , B2) ) 


> 


E3 


=multiply (Pro (z3) , -multiply (Pro (Pi/2) , B2) +multiply (Pro (0) , A2) ) 


> 


E4 


=multiply (Pro (z4) ,B1) : 


> 


EE 


=sum(sum(sum(sum( (El [ i] *E2 [ j ] *E3 [k] *E4 [1] ) , i=l . . 2) , j=l . . 2) , k=l 


> 


TrJ:=seq(sum(sum(op(v,EE) * (2-e) ,n=0. .1) ,m=0. .1) " (1+e) , v=l. . 16) : 


> 


zl 


=c*Pi/2+a*t: z2 : =d*Pi/2+b*t : z3 :=d*Pi/2+t : z4 :=c*Pi/2+t : 


> 


In 


=sum(TrJ[q] ,q=l . . 16) : 


> 


plot(ln(t),t=0. .Pi) : 



FIG. 3: A Maple implementation of Eq. (|l|). The polarizer regimes are encoded by setting the values of a,b,c and d to ±1, 0. 
't' is the skew angle. Setting the value of e to 1, adds the EM field contribution from the various centers before 'squaring;' 
as physics, this accounts for "crossover" signals that arrive simultaneously with the reflected signals and interfere. Setting e 
to zero, on the other hand, prevents interference of these signals. (This works with this model because of the peculiarities of 
down-conversion and a PBS; care must be taken modeling other setups to (in-)exclude appropriate signals before squaring.) 



tons,' or in classical terms: the signals, at their origin 
and simply carried along thereafter. 

The actual problem with entanglement arises elsewhere 
— with particle beams. A wave function describing 
a particle beam can not be considered simply as per- 
taining to a physical ensemble, because particles, one 
by one, are diffracted at slits, for example; they suf- 
fer, seemingly, 'entanglement' between Gibbsian ensem- 
ble members. Rationalizing this phenomenon requires 
other arguments^]; but, the desideratum of uniformity 
of principles, has lead to the mandate that radiation, too, 
be considered ontologically ambiguous until the moment 
of measurement, even though there is no need to do so. 
Its apparent digitization can be seen simply as a mani- 
festation of the nature of photoelectron detectors. 



Of course, the tactic employed herein works for all phe- 
nomena for which there is no Heisenberg Uncertainty, 
e.g., simple EPR correlations. |0 It offers a decidedly less 
mystical interpretation of many phenomena typically as- 
cribed to QM. Teleportation, for example, admits a pas- 
sive interpretation involving no "portation" of any nature. 
In the above model, so-called teleported states (1 and 4) 
are those which, although from separate random sources, 
eventually match up, and this is signaled by an appropri- 
ate coincidence between each's partner (2 and 3). Such 
an explanation is decidedly less enchanting than that con- 
veyed by the term "teleportation," but hugely more re- 
spectful of principles desirable for a rational explanation 
of the natural world. 
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